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Inelastic neutron scattering in variable magnetic field and high-field magnetization measurements, 
at the milikelvin temperature range, were performed to gain insight into the low-energy magnetic 
excitation spectrum and the field-induced level crossings in the molecular spin cluster {Crs}-cubane. 
These complementary techniques provide consistent estimates of the lowest level-crossing field. The 
overall features of the experimental data are explained using an isotropic Heisenberg model, based 
on three distinct exchange interactions linking the eight Cr m paramagnetic centers (spins s = 3/2), 
that is supplemented with a relatively large molecular magnetic anisotropy term for the lowest S = 1 
multiplct. It is noted that the existence of the anisotropy is clearly evident from the magnetic field 
dependence of the excitations in the INS measurements, while the magnetization measurements are 
not sensitive to its effects. 

PACS numbers: 75.25,+z, 75.50.Ee, 75.75.+a, 78.70.Nx 



I. INTRODUCTION 

Magnetic molecules are novel, physically realizable 
systems for exploring magnetic phenomena in low- 
dimensional magnetic materials^— . Although these sys- 
tems are very diverse, they possess the common charac- 
teristic of being achieved as crystalline samples of identi- 
cal molecules, each containing a relatively small number 
of mutually interacting paramagnetic centers ("spins"). 
In most cases, intermolecular magnetic interactions are 
negligible, as compared to intramolecular exchange inter- 
actions, so that experiments on bulk samples primarily 
probe the properties of the common, individual molecular 
unit. The small number of interacting spins in the basic 
molecular unit would suggest that it should be possible 
to formulate a theoretical model that can successfully ra- 
tionalize experimental results. However, as illustrated in 
the present work, in practice this task can present a sig- 
nificant challenge, especially when the experimental data 
are provided by complementary techniques. 

In the present study we explore a {Crg}-cubane mag- 
netic molecule, featuring eight interacting Cr m param- 
agnetic spin centers (spins s — 3/2)£~— . This system was 
also used to form a film on solid support by transfer of 
a monolayer at the air/ water interface by the Langmuir- 
Blodgett technique^. We report here detailed exper- 
imental results for several of the low-lying energy lev- 
els of this system as obtained by two complementary 
techniques, low temperature measurements of the mag- 



netic field dependence of the magnetization and inelas- 
tic neutron scattering (INS). We find that the magne- 
tization data can be accurately reproduced by a model 
based on purely isotropic Heisenberg exchange interac- 
tions between the eight Cr spins, but the field-dependent 
INS data are clearly inconsistent with a purely isotropic 
model. Previous INS studies have demonstrated the ef- 
fectiveness of the INS technique in directly determining 
the magnetic excitation spectrum of various magnetic 
molecules^—. Hence we introduce anisotropy effects 
in a simplified manner, which is successful in explain- 
ing these measurements, while maintaining results that 
are consistent with the magnetization measurements. 



In Sec. |TT1 various experimental details of the present 
work are described. These include a summary of the 
chemical synthesis of the {Cr§}-cubane compound in 
Sec. Ill Al and details of the experimental techniques are 
given for the magnetization and INS measurements in 
Sees. Ill Bl and III CI respectively. Our experimental re- 
sults are presented in Sec. IIIII and they are analyzed in 
terms of two theoretical models in Sec. IIVI Specifically, 
the magnetization results of Sec. IIII Al give rise to the 
isotropic Heisenberg model that is presented in Sec. II V Al 
and the INS results of Sec. IIII Bl lead us to supplement 
that model as described in Sec. IIVBI Finally, in Sec. |V] 
we summarize our present results and discuss open issues. 
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pies were placed in the core of an inductive coil of a self- 
resonant LC tank circuit, powered by a tunnel diode. The 
shift, Af(H), in the resonance frequency [11] of the coil 
due to the sample is proportional to the differential sus- 
ceptibility dM/dH. Measurements of A/ versus H were 
performed at fixed temperatures down to 300 mK using 
the 3 He cooling system with 18 Tesla superconducting 
magnet at the National High Magnetic Field Laboratory 
(NHMFL). The sweep rate was 0.3 T/min., and no hys- 
teresis effects were observed between sweeping up and 
down the magnetic field. 



FIG. 1. (color online) (a) Structure of the {Crs}-cubane mag- 
netic cluster. Phenyl groups and terminal benzoate ligands 
have been omitted for clarity, (b) scheme of the connectivity 
of the eight spin centers via superexchange pathways used in 
Sec. ITVAl 

II. EXPERIMENTAL DETAILS 

A. Chemical Synthesis and Characterization 

The {Cr 8 }-cubane cluster [Cr 8 04(02CPh)i 6 ] and its 
dcutcratcd analogue were synthesized following a modi- 
fied protocol of a previously published procedure^ em- 
ploying a crystallization method that uses fluoroben- 
zene and acetonitrile as solvents in order to reli- 
ably obtain high-quality crystalline material 9 -. Sam- 
ples were characterized by single-crystal X-ray diffrac- 
tion and IR spectroscopy. We note that there are min- 
imal structural differences between the present {Crs}- 
cubane compound [Cr 8 04(0 2 CPh)i 6 ]-4CH 3 CN-2H 2 
(space group C2/c) and the previously published 
analogue 8 [Crs04(02CPh)i6]-2Crl2Cl2 (space group 
P2\/c) due to crystal packing effects. 



B. Magnetization Measurements 

Values of the magnetization M versus external mag- 
netic field H for several low temperatures were obtained 
using two different measurement techniques, one based 
on a standard inductive method, and the second us- 
ing a tunnel diode oscillator (TDO) method. Utilizing 
fast digitizers at the Tohoku University High Magnetic 
Field Laboratory, the inductive method provided data 
for dM/dt and dH/dt which were integrated to give re- 
sults for M versus H for asymmetric half-cycle sweeps 
of duration 10 ms. Samples were immersed in liquid 
4 He and 3 He to maintain good thermal contact with the 
thermal bath. The resulting data for M versus H as ob- 
tained for the up and down portions of the half cycle were 
in good agreement, indicating that hysteresis effects are 
negligible. The measurements were performed on sam- 
ples obtained from four different synthesis batches, and 
for several temperatures in the range 0.42 to 1.5 K. For 
the TDO measurements polycrystalline protonated sam- 



C. Neutron Scattering Measurements 

Our INS measurements were performed both on 
deuterated and non-deuterated samples, with a base 
temperature of 60 mK. The experiments reported here 
were performed on approximately 2.6 g of polycrystalline 
deuterated samples to reduce attenuation and incoherent 
scattering from hydrogen atoms. Samples were sealed in 
a copper holder under helium atmosphere. The inelas- 
tic neutron spectra were collected using three different 
time-of-flight spectrometers: DCS at NCNRi& OSIRIS 
at ISIS, 19 and BASIS at SNS, ORNL3 For the mea- 
surements at the DCS we used a fixed incident energy 
Ej = 2.272 meV that yielded an energy resolution at 
zero energy transfer with a full width at half-maximum 
(FWHM) of approximately 60 /ieV. Additional data were 
collected at the OSIRIS back-scattering spectrometer, 
with a fixed final neutron energy E/ = 1.845 meV se- 
lected by a pyrolytic graphite PG(002) analyzer. The en- 
ergy resolution that this instrument provided at zero en- 
ergy transfer was approximately 20 /ieV. The third set of 
high-resolution inelastic neutron measurements were per- 
formed using the BASIS spectrometer that uses backscat- 
tering neutron reflections from Si(lll) analyzer crystals 
to select the final energy of the neutron of 2.08 meV. This 
instrument provides a ~ 4 fieV energy resolution at the 
elastic position. The BASIS data treatment involved sub- 
traction of a background spectrum measured separately 
with an empty Cu holder. The INS results reported here 
are all from OSIRIS, as they are the more detailed ones, 
however the partial results from other instruments are 
consistent with the findings on OSIRIS. 

III. EXPERIMENTAL RESULTS 

A. Level crossings from magnetization 

Shown in Fig. [2ja) is the resonance frequency of the 
tunnel diode oscillator for magnetic fields up to 18 T as 
measured at 300 mK. The high sensitivity of the TDO 
method is demonstrated by the observation (inset of Fig. 
[2ja)) of the sharp peak at approximately 4.8 T that can 
be identified with the proton magnetic resonance absorp- 
tion line. Shown in Fig. &h) is the field dependence of 
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the frequency shift A/, equivalently dM/dH, due to the 
intrinsic Cr 8 -cubane sample, as obtained by subtracting 
the background signal seen in Fig. 3(a) associated with 
an empty resonator. 

Level crossing fields were identified by the appearance 
of peaks in the differential susceptibility dM/dH versus 
H as measured at fixed low temperatures, and these data 
are shown in Table Q] Specifically, columns 1-4 of the 
Table list the values of the level-crossing fields that were 
obtained using the inductive method. (These data were 
collected using four different samples, with fields in the 
range from to 20 T.) The corresponding field values as 
obtained using the TDO method up to 18 T for a fifth 
sample are listed in column 5. The TDO peak at 17.9 
T is only partially covered by the field sweep. It should 
be emphasized that the two techniques provide similar 
results. 
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FIG. 2. Variation with the magnetic field of the resonance 
frequency in a tunnel-diode oscillator, (a) Raw data measured 
at T = 0.3A". Inset shows the sharp line due to proton NMR 
at approximately 4.8 T. (b) The frequency shift of the intrinsic 
signal at T = 0.3 K after background subtraction. 



B. Inelastic neutron scattering 

Figure [3] shows intensities integrated over the Q = 
0.3 to 1.8 A -1 range vs neutron energy loss, on a 
semilog plot at the lowest temperature achievable with 
the cryo-magnet, nominally T = 60 mK. The large 
elastic peak has a FWHM of ^0.021 meV consistent 



TABLE I. Comparison of the measured level crossings fields 
for five different samples and the prediction of the isotropic 
Heisenberg model (see Eq. Q). Columns 1-4: pulsed- field 
method; column 5: TDO. The estimated error for all the 
pulsed-field entries is 0.15 T. For the TDO data the errors, 
listed in parenthesis, were obtained by a Lorentzian-based 
multi-peak fitting procedure. 
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FIG. 3. (color online) (a) Inelastic neutron scattering spec- 
trum at a nominal temperature of 60 mK. The intensity axis 
uses a logarithmic scale. The dashed line shows a best esti- 
mate of the nonmagnetic background (see text). Inset: The 
background subtracted scattering consists of three excitations 
at approximately 0.16, 0.32, 0.48 meV with similar widths, 
and a much broader and weaker excitation centered at ~0.85 
meV. (The error bars in all figures represent one standard 
deviation in the measurement.) 



with the expected instrumental energy resolution for this 
configuration^ Two broader peaks are clearly visible 
as superpositions on the otherwise bell-shaped resolu- 
tion function. To obtain the genuine magnetic spectra 
from the sample we reduce the data using a procedure 
we employed in a similar INS study of the {Mo72Fe3o} 
magnetic molecule^ The quasi-elastic term representing 
the instrumental resolution function, the incoherent and 
static disorder due to the sample (we name background 
curve), is described by a sum of two co-centered peaks: 
a dominant Gaussian and a minor Lorentzian. To ob- 
tain reliable parameters of the Gaussian/Lorentzian, the 
intensities at various temperatures are refined simulta- 
neously, with extra peaks due to the magnetic spectra, 
while maintaining the same values of these parameters for 
the refinements of the various data sets at all measured 



temperatures. The resulting background curve (that in- 
cludes non-magnetic contributions from the sample and 
holder) is shown as a dashed line in Fig. |3l Subtracting 
the background curve yields the genuine magnetic spec- 
tra as shown in the inset of Fig. [3j The detailed analysis 
shows that the spectra at base temperature consists of 
three excitations at approximately 0.16, 0.32, 0.48 meV 
with similar widths, and a much broader and weaker ex- 
citation centered at ^0.85 meV. 
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FIG. 5. (color online) Magnetic spectra at T = 0.06K at 
various applied magnetic fields obtained from INS measure- 
ments (shifted for clarity). The spectra are obtained after the 
subtraction of the background that includes the quasi-elastic 
non-magnetic contributions from the sample and the instru- 
mental resolution, in a procedure similar to the one employed 
to obtain the inset of Fig. [3] The dashed lines show the pro- 
gression of excitations with the increase of magnetic field. 



FIG. 4. (color online) Magnetic excitations at various tem- 
peratures as indicated (shifted for clarity). The data are ob- 
tained after background subtraction that includes the quasi- 
elastic non-magnetic contributions from the sample and the 
instrumental resolution, using a procedure similar to the one 
employed to obtain the inset of Fig. [3] 

Figure 0] displays the magnetic excitation spectra at 
various temperatures after the subtraction of the back- 
ground curve as described above. The general effect of 
increasing temperature is in reducing the intensities of 
the excitations. However, above sa 3 K, the lowest ex- 
citation broadens and shifts to lower energies, whereas 
the other excitations broaden and almost disappear. The 
broadening and shift to lower energy of the lowest excita- 
tion as the temperature is raised is commonly understood 
to be due to the occupation of excited states with access 
to many more allowed transitions some of which are asso- 
ciated with smaller gaps than the first excited states. It 
should be noted that the experimental spectra at 60 and 
800 mK are, within uncertainty, identical. We therefore 
argue that a comparison of the INS measurements at 60 
mK with the magnetization measurements at T ~ 0.3 K 
is justified. 

To confirm the magnetic origin of the excitations and 
to get more insightful characteristics of the spectra, 
we conducted the INS experiments under applied mag- 
netic field at various temperatures. The background- 
subtracted spectra at T = 0.06 K are shown in Fig. [S] 



(a) and (b). Similar results are obtained at T = 0.8 
K. The lowest excitation (0.16 meV) is the most sensi- 
tive to the magnetic field. This excitation shifts to lower 
energies and disappears between 2 and 3 T, and then 
reappears above ~ 3 T. This can be interpreted as a di- 
rect observation of level crossing at approximately 2.8 T. 
This is consistent with the first level crossing observed in 
the magnetization measurements (» 2.8 T). In fact, INS 
raw data collected on BASIS (SNS) at various applied 
magnetic fields shows the same level crossing of the low- 
est excitation. The high resolution setup (4 /LteV) allows 
one to clearly observe the level crossing, indicated by the 
reopening of the energy gap (Fig. [5]). 

The behavior of the lowest excitations as a function of 
magnetic field is shown in Fig. |Bfb). The excitation at 
0.32 meV practically remains unchanged and broadens 
in the vicinity of the level crossing, and shifts to higher 
energies with the increase of magnetic field when the new 
ground state is established following the level crossing 
(traced as guide in Fig. OJb)). The 0.48 meV excitation 
shifts to higher energy with increasing magnetic field and 
beyond the level crossing increases at a faster rate (shown 
in dashed line as guide in Fig.[5jb)). At higher fields (> 4 
T), excitations that are not accounted for by the S = 0, 
and S = 1 energy levels emerge. We hypothesize that 
they arise from excitations that involve the lowest 5 = 2 
multiplet. 
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FIG. 6. (color online) INS raw data collected on BASIS (SNS) 
at various applied magnetic fields at T — 1.5 K. The high 
resolution setup (4 /ieV) allows direct observation of the level 
crossing with increasing field, as the excitation shifts towards 
zero energy (below 2 Tesla) and the energy gap reopens at 
higher fields. 



IV. COMPARISON WITH THEORY 



A. Isotropic Model 



As a first attempt towards building a comprehensive 
theoretical model yielding results consistent with our 
experimental data, we utilized an isotropic Heisenberg 
Hamiltonian based on an initial assumption that the 
Cr m ions in their octahedral CrOg ligand field envi- 
ronments should be well approximated by pure spin-3/2 
centers 2 ^. Three distinct exchange constants are used 
to describe the interactions between the eight Cr m spin 
centers as shown in Fig. lb. (Contrary to an earlier 
attemptiS in which pathways involving one /i-oxo and one 
/i-carboxylate group were distinguished from pathways 
involving one /i-oxo and two /i-carboxylate groups, we 
here operate with a simplified model in which these ex- 
change pathways are characterized by the same exchange 
constant.) The Hamiltonian, %, of the spin system in the 



presence of an external magnetic field H is given by 

% =J\ (Si • S 2 + S; 2 • S. 3 + S*2 • S*4 + S. 3 • S i 

+ £3 -£4 + £4 -£i) 
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+ J 3 (si • 5*7 + S4 ' S 5 + S 2 ' S 8 + S3 ' S S ) 
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where the individual spin operators are those for intrin- 
sic spins s — 3/2, they are given in units of h, and we 
used the value g — 1.985 for the spectroscopic splitting 
factor appropriate for Cr m spin centers. (Tildes, written 
below symbols, are used to denote quantum operators.) 
Inasmuch as the total spin operators 5 2 and S z commute 
with H, the eigenstates of these operators are described 
by quantum numbers 5 and M$ whose values range from 
to 12 and from —5 to 5, respectively. Each multiplet 
has (25 + l)-fold degeneracy when H = 0. In an exter- 
nal magnetic field the degeneracy is lifted due to a shift, 
g^sH Ms, originating from the Zeeman term of Eq. (fl"j). 
As the external field is increased from H = 0, the ground 
state changes successively from 5 = 0, Ms = to 5 = 1, 
A/5 = —1, etc., in integer steps of 5 and Mg. Each 
of these changes of the ground-state quantum numbers 
is referred to as a ground state level crossing, and the 
field at which the ground state changes from 5 — 1 to 
5 will be denoted by Hs . It is easy to show that the 
values of the Hs can be determined using the relation 
Es — Es-i = gusHs, (5 = 1, 2, 3, ...), where Es denotes 
the zero-field energy of the lowest-energy multiplet with 
quantum number 5. (It is convenient to choose Eq =0.) 

To diagonalize the Hamiltonian of Eq. fl} we used the 
numerical package MAGPACK, 22 so as to provide the 
full set of energy levels as well as the temperature de- 
pendence of the weak-field susceptibility. The values of 
the exchange constants were first constrained using the 
method described in Ref . |_10| so as to reproduce the weak- 
field susceptibility (measured at 0.5 T using a Quantum 
Design MPMS magnetometer and shown in Fig.[7j) above 
100 K. We then varied the values of the constrained ex- 
change constants so as to optimize the fit between the 
calculated and measured values of the first five ground 
state level-crossing fields (see Table I). The resulting opti- 
mized values of the exchange constants are Ji/Ub = 32.1 
K, J 2 /k B = 2.2 K, and J 3 /k B = -23 K. The correspond- 
ing theoretical values of the ground state level-crossing 
fields are listed in the right-most column of Table HI and 
it is noted that they are indeed in good agreement with 
our measured values 

For the optimized isotropic Hamiltonian model the to- 
tal angular momentum quantum number of the ground 
state for zero magnetic field is 5 = 0. This result is 
consistent with what can be inferred from the experi- 
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FIG. 7. (color online) Molar magnetic susceptibility, \ — 
M/H, of {Crg}-cubane for 0.5 T: Experiment (open cir- 
cles) and theory (solid curve) as obtained using the isotropic 
Heisenberg model, Eq. fT]) for optimized values of the ex- 
change constants. The approach of T\ — > as T — > shown 
in the inset indicates that the ground state of the system is 
S — (see text for more detail). 



mental weak-field susceptibility data. We recall the zero- 
field fluctuation formula that the limiting value of T\ for 
T -» is given by N A g 2 ^ 2 B S{S + l)/(3fc B ) for a ground 
state with quantum number S. (Here Na is Avogadro's 
number and kg is Boltzmann's constant.) Indeed, the 
data for T%, shown in the inset of Fig. tend towards 
zero for T — > 0, whereas if S = 1 for the ground state 
the limiting low-temperature value would be T\ m 1 cm 3 
K/mol. 

For completeness, we list here the values of the energy 
(in meV) and the ^-quantum number for those levels of 
the isotropic model lying within 2 meV of the ground 
state: 0.338 (S = 1), 1.025 (S = 2), 1.173 (S = 1), 1.291 
(S = 0), 1.660 (S = 1), 1.916 (S = 2) meV. 

Finally, we remark that the susceptibility that is cal- 
culated using the optimized isotropic Hamiltonian (solid 
curve in Fig. [7J does not provide a good fit to the mea- 
sured weak-field susceptibility at low temperatures. This 
fact and the field dependence of the INS data strongly 
suggest that non-Heisenberg terms play a significant role 
in the description of {Crs}-cubane. 



B. Including zero-field splitting 

Although the isotropic Hamiltonian of Eq. (JTJ) is able 
to accurately reproduce the measured ground state level- 
crossing fields, our INS measurements provide clear evi- 
dence of additional non-Heisenberg interactions. Namely, 
as remarked above, the energy spectrum obtained from 
Eq. ([1]) will have a lowest S = 1 multiplet that is three- 



fold degenerate (with energy E\) for H = 0, that splits 
into three levels {E-i,E\ ± g^sH) for H ^ 0. Such be- 
havior would be detectable in an INS measurement as a 
single peak for H = 0, that would split into three satel- 
lite peaks for H > 0. However, as seen in Fig.O our INS 
data show the presence of three (not one) distinct peaks 
in zero field that, for increasing field strengths, consist 
of 1) a low energy satellite, starting from 0.16 meV, that 
shifts to lower energy, 2) a high energy satellite, start- 
ing from 0.48 meV, that shifts to higher energy, and 3) a 
central peak at 0.32 meV. 

We show now that the observed behavior of these three 
INS peaks in a magnetic field can be qualitatively un- 
derstood in terms of an S — 1 triplet whose degener- 
acy in zero field has been lifted by an anisotropy mecha- 
nism rooted in the coordination environment of the four 
outer CrC>6 octahedra. Specifically, a closer inspection of 
the individual chromium ligand field environments in the 
molecular structure of Crg-cubane reveals that two types 
of coordination environments need to be differentiated 9 . 
Although the environment of the central four CrC>6 oc- 
tahedra shows only relatively small deviations from the 
ideal O/j symmetry (Cr-0 bond distances only vary by 
2%, O-Cr-0 bond angles by less than 10° from the 90° 
ideal), the outer CrC>6 groups experience a stronger dis- 
tortion, primarily caused by terminal benzoate ligands, 
leading to very small (65°) and very large (108°) O-Cr- 
O bond angles. This significant distortion in turn re- 
sults in single-ion zero-field splitting that we quantify us- 
ing a phenomenological anisotropy model. Notably, the 
molecular symmetry clearly deviates from the tetrahe- 
dral ideal: The crystallographically imposed point group 
symmetry of the Crs-cubane molecules is Ci . The sub- 
structure of the eight Cr centers for example differs from 
a fully Trf-symmetric arrangement by more than 0.48 A. 
In this absence of tetrahedral or cubic symmetry, the 
single-ion anisotropy tensors of the Cr spin centers thus 
add up to a non-zero molecular anisotropy tensor, which 
supports the principal assumption of the molecular zero- 
field splitting. We thus supplement the isotropic Hamil- 
tonian of Eq. (fT]) with a general molecular anisotropy 
term of the form 

%' = E 1 l + g^ B H-S + d{S 2 z - h + e(S 2 x - Si), (2) 

3 u 

which we apply to the reduced (three-dimensional) sub- 
space of the lowest S = 1 triplet. The first term in Eq. ([2]) 
denotes the product of the field-free multiplet energy E\ 
with the (3 x 3) unit matrix, 1_; the Cartesian axes denote 
symmetry axes of a given molecule; and the three Carte- 
sian components of the vector operator S are given by the 
standard (3x3) matrix representation of angular momen- 
tum operators for the S = 1 manifold. For zero magnetic 
field the eigenvalues of H' are E\ — 2rf/3, E\ + d/3 ± e. 
We adopt the values E x = 0.32 meV, d = 3e = 0.24 meV 
thereby reproducing the energies (0.16 meV, 0.32 meV, 
and 0.48 meV) of the zero-field peaks that are observed 
in our INS measurements. 
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FIG. 8. (color online) (a) Peak INS energies (symbols) and 
excitation bands (shaded) versus magnetic field. The exci- 
tation bands were obtained from Eq. ((2J as described in the 
text, (b) Energy bands versus magnetic field, obtained from 
Eq. as described in the text. 



Although the three zero-field eigenvalues of %' are in- 
dependent of a molecule's orientation, this is not the case 
for H > 0. Therefore, in order to describe an actual 
polycrystalline sample of {Crsj-cubane in the presence 
of a magnetic field, we have diagonalized %' for many 
(10 6 ) uniformly distributed orientations of the Cartesian 
axes for each value of H. Each orientation produces a 
different set of three eigenvalues, and as the number of 
orientations becomes very large, these form three con- 
tinuous distributions (bands) of energies. To show the 
variation of these distributions with field, we have calcu- 
lated the mean value of each distribution plus (minus) its 
standard deviation, which define the energy bands that 
are shown in Fig.^b). It can be seen that the lowest en- 
ergy band intersects the S = ground state in the range 
2.2 T < H < 2.8 T; i.e., for some orientations the level- 
crossing field is as small as 2.2 T, and for other orienta- 
tions it is as large as 2.8 T. This range of fields compares 
favorably with the 2.6 - 3.1 T range of observed values 
of the first level-crossing field, that are listed in Table I. 
Since Eq. ([2]) does not produce a unique level-crossing 
field, this suggests that at very low temperatures peaks 
in dM/dH should be broader than the peaks that are ob- 
served in a system without appreciable anisotropy, while 



at higher temperatures, this effect will be hidden within 
the thermal broadening. The width of the peaks shown 
in Fig. are consistent with the width as calculated 
using the isotropic model, Eq. (JXJ) , for T — 300 mK; but 
we expect that the measured widths would reveal the ef- 
fect of anisotropy if the experiment were repeated at a 
much lower temperature. 

As described above, we have used Eq. ^ to obtain 
three field-dependent energy eigenvalues for each of 10 6 
different molecular orientations. Including the (field- 
independent) 5 = eigenvalue, this gives 4 x 10 6 eigen- 
values for each value of magnetic field strength. In order 
to compare these results with the field-dependence of the 
INS measurements, we note that at 60 mK (the temper- 
ature used for the measurements shown in Fig. [5]) only 
the ground state of a given molecule is expected to have 
significant thermal occupation, so transitions between ex- 
cited states should be negligible. Therefore, we have cal- 
culated the energy differences between the ground state 
and the three excited states for each of the 10 6 orienta- 
tions and for each value of magnetic field strength. (As 
described above, the ground state is 5 = for H < 2 
T, it is 5 = 1 for H > 3 T, and the value of the level 
crossing field is different for different orientations.) Be- 
cause these energy differences depend on the orientation 
of a molecule, they generate three excitation bands, and 
these are shown in Fig.JSJa), along with the field depen- 
dence of the INS peaks. It can be seen that all three 
excitation bands exhibit the same qualitative behavior 
as observed for the INS peaks. Specifically, the lowest 
band decreases toward zero for H < 2 T, then increases 
for H > 3 T; the middle band is roughly independent 
of field for H < 3 T, then slowly increases with field for 
H > 3 T; and the highest band has a positive slope for all 
H > 0. The poorer agreement visible above 4 T may be 
due to zero-field splitting effects of the lowest 5 = 2 mul- 
tiplet. These effects become significant for field strengths 
in the vicinity of the second level crossing field (« 6 T). 

Finally, we have also considered how the temperature 
dependent magnetic susceptibility is affected by the ob- 
served zero-field splitting. Our data for the weak-field 
(0.5 T) susceptibility (Fig. 7), extending down to 2 
K, appear to be reaching a maximum (0.28 cm 3 /mol) 
somewhere in the temperature interval (0, 2) K. By con- 
trast the theoretical susceptibility as obtained using the 
optimized isotropic Heisenberg model has a maximum, 
0.25 cm 3 /mol, that occurs at the much high temperature 
T = 11.6 K. As remarked earlier, we take this as the 
second piece of evidence, besides our INS findings, for 
the basic inadequacy of a strictly isotropic description 
of the system. To estimate how the anisotropy terms 
of Eq. (2) influence the low-temperature (< 5 K) sus- 
ceptibility, we have calculated M/H versus T using Eq. 
(2) including only the four lowest energy levels that are 
discussed above (energies 0, 0.16, 0.32, 0.48 meV). The 
effect of anisotropy is very large: The resulting suscepti- 
bility has a maximum, 0.16 cm 3 /mol, at 2.1 K. While it is 
gratifying that the peak temperature is consistent with 
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that of the available experimental data, the low value 
of the maximum is unsatisfactory. We have attempted 
to resolve this discrepancy by assuming the presence of 
independent paramagnetic impurities, i.e., by adding a 
Curie term, C/T, to the susceptibility of the four-level 
model. (Specifically, we chose C = 0.25 cm 3 -K/mol so 
as to achieve agreement between the theoretical and ex- 
perimental values (0.28 cm 3 /mol) of the susceptibility for 
T = 2 K.) This attempt was unsuccessful since the re- 
sulting theoretical value for T — 5 K is 0.17 cm 3 /mol, 
as compared to the experimental value 0.27 cm 3 /mol. 
The inadequacy of the four-level model may be due to 
not having included the other low-lying energy levels, in 
particular the anisotropy-perturbed levels of the lowest 
5 = 2 multiplct. Unfortunately, at the present time the 
form and the strength of those anisotropy terms are un- 
known. 



V. SUMMARY 

The present article reports an analysis of the magnetic 
field dependence of the energy spectrum for the {Cr 8 }- 
cubane magnetic molecule using multiple experimental 
techniques. Ground state level crossings were observed 
using both tunnel diode oscillator (TDO) measurements 
and pulsed-field magnetization measurements, and ex- 
citation energies were measured directly using inelastic 
neutron scattering (INS) spectroscopy. In particular, 
peaks in dM/dH were observed in both the TDO and 
pulsed-field measurements, which we identify as ground 
state level-crossing fields. These field values (given in 
Table I) are used in conjunction with high temperature 
susceptibility data to determine the exchange constants 
for the isotropic Heisenberg Hamiltonian of Eq. ([1} . 

We note that the TDO and pulsed-field measurements 
provide no hint of the need for anything other than an 
isotropic Heisenberg model, while it is clear from our INS 
measurements that such a model is incomplete. Specif- 



ically, we observed three INS peaks, and based on their 
field dependence, we are able to identify the energy lev- 
els as belonging to the lowest 5 = 1 multiplet. These 
peaks show that this lowest S = 1 multiplet is split, even 
in zero magnetic field, which is unexplainable in terms 
of the Heisenberg Hamiltonian.— The magnitude of the 
zero-field splitting (0.16 meV) is comparable to that ob- 
served for the antiferromagnetic ring [CrsFsPivig]^ 

We attribute the observed field dependence of the INS 
data as being due to the effects of anisotropy not in- 
cluded in the initial (Heisenberg) Hamiltonian, and we 
hence supplement that Hamiltonian with an additional 
perturbing term, given in Eq. @, that incorporates the 
zero-field splitting of the lowest 5 = 1 multiplet. We 
find that by including this anisotropy term, we are able 
to achieve a semi-quantitative understanding of the INS 
data, while maintaining results that are consistent with 
the magnetization measurements. 
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